Chapter 9 

Post-main sequence evolution through 
helium burning 


To describe the evolution of stars after the main sequence, it is useful to make a division based on the 
mass: 

low-mass stars are those that develop a degenerate helium core after the main sequence, leading to a 
relatively long-lived red giant branch phase and to the ignition of He in a so-called helium flash. 
This occurs for masses between 0.8 M© and ~ 2 M© (this upper limit is sometimes denoted as 

^HeF)- 

intermediate-mass stars ignite helium stably in a non-degenerate core, but after central He burning 
they develop a degenerate carbon-oxygen core. They have masses between Af Hc F and M up « 
8 M q . Both low-mass and intermediate-mass stars shed their envelopes by a strong stellar wind 
at the end of their evolution and their remnants are CO white dwarfs. 

massive stars have masses larger than M up » 8 M© and ignite carbon in a non-degenerate core. 
Except for a small mass range (« 8 - 11 Af©) these stars also ignite heavier elements in the core 
until an Fe core is formed which collapses. 


9.1 The Schonberg-Chandrasekhar limit 


During central hydrogen burning on the main sequence, we have seen that stars are in thermal equi¬ 
librium (r nuc » trh) with the surface luminosity balanced by the nuclear power generated in the 
centre. After the main sequence a hydrogen-exhausted core is formed inside which nuclear energy 
production has ceased. This inert helium core is surrounded by a hydrogen-burning shell and a H-rich 
envelope. For such an inert core to be in thermal equilibrium requires l(m ) = J m e nuc dm - 0 and 
hence dT/dr oc / = 0, which implies that the core must be isothermal to remain in TE. Such a stable 
situation is possible only under certain circumstances. 

If ideal-gas conditions hold in the core, then an upper limit exists for the mass of the core M c 
relative to the total mass of the star. If q c - MJM exceeds this limit, then the pressure within an 
isothermal core cannot sustain the weight of the overlying envelope. This was first pointed out by 
Schonberg and Chandrasekhar in 1942, who computed the limiting core mass fraction to be 


9< 9sc = 0.37(^)- 

M v li c ' 
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(9.1) 




where ju c and /z env are the mean molecular weight in the core and in the envelope respectively. This 
limit is known as the Schonberg-Chandrasekhar limit. The typical value qsc ~ 0.10 is appropriate 
for a helium core with /r c = 1.3 and a H-rich envelope. (A simple, qualitative derivation of eq. 9.1 
can be found in K&W section 30.5.1.) 

Stars that leave the main sequence with a helium core mass below the Schonberg-Chandrasekhar 
limit can therefore remain in complete equilibrium (HE and TE) during hydrogen-shell burning. This 
is the case for stars with masses up to about 8 M 0 , if convective overshooting is neglected. Over¬ 
shooting increases the core mass at the end of central H-buming, and therefore the upper mass limi t 
for stars remaining in TE after the main sequence decreases to about 2 M e in calculations that include 
moderate overshooting. 

When the mass of the H-exhausted core exceeds the Schonberg-Chandrasekhar limit - either im¬ 
mediately after the main sequence in relatively massive stars, or in lower-mass stars after a period 
of H-shell burning during which the helium core mass increases steadily - thermal equilibrium is no 
longer possible. The helium core then contracts and builds up a temperature gradient. This tempera¬ 
ture gradient adds to the pressure gradient that is needed to balance gravity and keep the star in HE. 
However, the temperature gradient also causes an outward heat flow from the core, such that it keeps 
contracting and heating up in the process (by virtue of the virial theorem). This contraction occurs on 
the thermal (Kelvin-Helmholtz) timescale in a quasi-static way, always maintaining a state very close 
to HE. 




Figure 9.1. Evolution tracks for stars of quasi-solar composition (X = 0.7, Z = 0.02) and masses of 1, 2, 3, 5, 
7 and 10M o in the H-R diagram (left panel) and in the central temperature versus density plane (right panel). 
Dotted lines in both diagrams show the ZAMS, while the dashed lines in the right-hand diagram show the 
borderlines between equation-of-state regions (as in Fig. 3.4). The 1 M e model is characteristic of low-mass 
stars: the central core becomes degenerate soon after leaving the main sequence and helium is ignited in an 
unstable flash at the top of the red giant branch. When the degeneracy is eventually lifted. He burning becomes 
stable and the star moves to the zero-age horizontal branch in the HRD, at logL « 1.8. The 2 M e model is 
a borderline case that just undergoes a He flash. The He flash itself is not computed in these models, hence 
a gap appears in the tracks. The 5 M e model is representative of intermediate-mass stars, undergoing quiet 
He ignition and He burning in a loop in the HRD. The appearance of the 7 and 10 M 0 models in the HRD 
is qualitatively similar. However, at the end of its evolution the 10 M e star undergoes carbon burning in the 
centre, while the cores of lower-mass stars become strongly degenerate. (Compare to Fig. 7.4.) 
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Low-mass stars (M < 2 M 0 ) have another way of maintaining both HE and TE during hydrogen- 
shell burning. In such stars the helium core is relatively dense and cool and electron degeneracy can 
become important in the core after the main sequence. In that case the Schonberg-Chandrasekhar 
limit no longer applies, because degeneracy pressure can support the weight of the envelope even in 
a relatively massive core. 1 Inside such degenerate helium cores efficient energy transport by electron 
conduction can keep the core almost isothermal. 

Effects of core contraction: the ‘mirror principle’ 

The following principle appears to be generally valid, and provides a way of interpreting the results 
of detailed numerical calculations: whenever a star has an active shell-burning source, the burning 
shell acts as a mirror between the core and the envelope: 

core contraction => envelope expansion 

core expansion => envelope contraction 

This result can be understood by the following argument. To maintain thermal equilibrium, the burn¬ 
ing shell must remain at approximately constant temperature due to the thermostatic action of nuclear 
burning. Contraction of the burning shell would entail heating, so the burning shell must also remain 
at roughly constant radius. As the core contracts, p s h e ii must therefore decrease and hence also the 
pressure in the burning shell must decrease. Therefore the pressure P e nv of the overlying envelope 
must decrease, so the layers above the shell must expand (an example of this behaviour can be seen 
in Fig. 9.3, to be discussed in the next section). 

9.2 The hydrogen-shell burning phase 

In this section we discuss in some detail the evolution of stars during hydrogen-shell burning, until 
the onset of helium burning. Based on the above section, qualitative differences are to be expected 
between low-mass stars (M < 2 M 0 ) on the one hand and intermediate- and high-mass stars (M > 
2 M q ) on the other hand. Therefore we discuss these two cases separately, starting with the evolution 
of higher-mass stars because it is relatively simple compared to low-mass stars. We use two detailed 
stellar evolution sequences, for stars of 5 M 0 and 1M 0 respectively, as examples for the general 
evolutionary behaviour of stars in these two mass ranges. 

9.2.1 Hydrogen-shell burning in intermediate-mass and massive stars 

Fig. 9.2 shows the evolution track of a 5 M 0 star of quasi-solar composition (A = 0.7, Z = 0.02) in the 
H-R diagram, and Fig. 9.4 shows some of the interior details of the evolution of this star as a function 
of time from the end of central hydrogen burning. Point B in both figures corresponds to the start of 
the overall contraction phase near the end of the main sequence (when the central H mass fraction 
X c ~ 0.03) and point C corresponds to the exhaustion of hydrogen in the centre and the disappearance 
of the convective core. The hatched regions in the ‘Kippenhahn diagram’ (lower panel of Fig. 9.4) 
show the rapid transition at point C from hydrogen burning in the centre to hydrogen burning in a 
shell. 

The H-exhausted core initially has a mass of about 0.4 M 0 which is below the Schonberg-Chandra- 
sekhar limit, so the star initially remains in TE and the first portion of the hydrogen-shell burning 

1 As long as the degenerate core mass does not exceed the Chandrasekhar mass - note the very different physical mean¬ 
ings of the Chandrasekhar mass and the Schonberg-Chandrasekhar limit. 
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Figure 9.2. Evolution track in the 
Hertzsprung-Russell diagram of a 5 M 0 star 
of initial composition X = 0.7, Z = 0.02. See 
text for details. 


log Teff(K) 


phase (C-D) is relatively slow, lasting about 2 x 10 6 yr. The temperature and density gradients be¬ 
tween core and envelope are still shallow, so that the burning shell initially occupies a rather large 
region in mass. This phase is therefore referred to as thick shell burning. The helium core gradually 
grows in mass until it exceeds the S-C limit and the contraction of the core speeds up. The envelope 
expands at the same time, exemplifying the ‘mirror principle’ discussed above. This becomes more 
clear in Fig. 9.3 which shows the radial variations of several mass shells inside the star. After point 
C the layers below the burning shell contract while the layers above expand, at an accelerating rate 
towards the end of phase C-D. As a result the temperature and density gradients between core and 
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Figure 9.3. Radial variation of various mass 
shells (solid lines) in the 5 M Q (Z = 0.02) 
star of Fig. 9.4, during the early post-main se¬ 
quence evolution. Each line is labelled with 
its mass coordinate m in units of M e ; the 
top-most curve indicates the total radius R. 
Gray areas indicate convection and red cross- 
hatched areas have intense nuclear burning 
(e nuc > 10 L/M). Letters B.. .E correspond 
to those in Fig. 9.4. 


age (10 6 yr) 
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Figure 9.4. Internal evolution of a 5 M e star 
of initial composition X = 0.7, Z = 0.02. The 
panels show various internal quantities as a 
function of time, from top to bottom: 

(a) Contributions to the luminosity from hy¬ 
drogen burning (red line), helium burning 
(blue) and gravitational energy release (or¬ 
ange; dashed parts show net absorption of 
gravitational energy). The black line is the 
surface luminosity. 

(b) Central mass fractions of various elements 
('H, 4 He, 12 C, 14 N and ie O) as indicated. 

(c) Internal structure as a function of mass 
coordinate m, known as a ‘Kippenhahn di¬ 
agram’. A vertical line through the graph 
corresponds to a model at a particular time. 
Gray areas are convective, lighter-gray ar¬ 
eas are semi-convective. The red hatched re¬ 
gions show areas of nuclear energy genera¬ 
tion, where e nu c > 10 L/M (dark red) and 
fnuc > 2 L/M (light red). The letters B.. .J 
indicate the corresponding points in the evo¬ 
lution track in the H-R diagram, plotted in 
Fig. 9.2. See text for details. 


age (10 6 yr) 


envelope increase, and the burning shell occupies less and less mass (Fig. 9.4c). The latter portion of 
hydrogen-shell burning is therefore referred to as thin shell burning. Most of the time between C and 
D is spent in the thick shell burning phase at relatively small radii and log T e g- > 4.05. The phase of 
expansion from log T e f[ ~ 4.05 to point D at log T e g & 3.7 occurs on the Kelvin-Helmholtz timescale 
and takes only a few times 10 5 yrs. A substantial fraction of the energy generated by shell burning is 
absorbed by the expanding envelope (dashed yellow line in Fig. 9.4a), resulting in a decrease of the 
surface luminosity between C and D. 

The rapid evolution on a thermal timescale across the H-R diagram from the end of the main 
sequence to r e g- « 5000 K is characteristic of all intermediate-mass stars. The probability of detecting 
stars during this short-lived phase is very small, resulting in a gap in the distribution of stars in the 
H-R diagram known as the Hertzsprung gap. 

As point D is approached the envelope temperature decreases and the opacity in the envelope rises, 
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impeding radiative energy transport. The envelope grows increasingly unstable to convection, starting 
from the surface, until at D a large fraction of the envelope mass has become convective. During 
phase D-E the star is a red giant with a deep convective envelope. The star is then located close to the 
Hayashi line in the H-R diagram, and while it continues to expand in response to core contraction, 
the luminosity increases as the effective temperature remains at the approximately constant value 
corresponding to the Hayashi line. The expansion of the star between D and E still occurs on the 
thermal timescale, so the H-shell burning phase of intermediate-mass stars on the red-giant branch is 
very short-lived. 

At its deepest extent at point E, the base of the convective envelope is located at mass coordinate 
m = 0.9 M© which is below the maximum extent of the former convective core during central H- 
buming (about 1.25 M© at the start of the main sequence). Hence material that was formerly inside 
the convective core, and has therefore been processed by hydrogen burning and the CNO-cycle, is 
mixed throughout the envelope and appears at the surface. This process is called dredge-up and 
occurs about halfway between D and E in Fig. 9.2. Dredge-up on the red giant branch also occurs in 
low-mass stars and we defer its discussion to Sec. 9.2.3. 

The helium cores of intermediate-mass stars remain non-degenerate during the entire H-shell 
burning phase C-E, as can be seen in Fig. 9.1. These stars develop helium cores with masses larger 
than 0.3 M©, the minimum mass for helium fusion discussed in Ch. 7. In the 5 M© star at point E the 
helium core mass is 0.6 M© when a central temperature of 10 8 K is reached and helium is ignited in 
the core. The ignition of helium halts further core contraction and envelope expansion and therefore 
corresponds to a local maximum in luminosity and radius. Evolution through helium burning will be 
discussed in Sec. 9.3.1. 

9.2.2 Hydrogen-shell burning in low-mass stars 

Compared to intermediate-mass stars, low-mass stars (with M <2 M©) have small or no convective 
cores during central hydrogen burning, and when they leave the main sequence their cores are rel¬ 
atively dense and already close to becoming degenerate (see Fig. 9.1). In stars with M < 1.1 M© 
the transition from central to shell hydrogen burning is gradual and initially MJM < 0.1 so the star 
can remain in thermal equilibrium with an isothermal helium core. By the time the helium core has 
grown to « 0.1 M, its density is large enough that electron degeneracy dominates the pressure and the 
Schonberg-Chandrasekhar limit has become irrelevant. Therefore low-mass stars can remain in HE 
and TE throughout hydrogen-shell burning and there is no Hertzsprung gap in the H-R diagram. 

This can be seen in Fig. 9.5 which shows the internal evolution of a 1 M© star with quasi-solar 
composition in a Kippenhahn diagram and the corresponding evolution track in the H-R diagram. Hy¬ 
drogen is practically exhausted in the centre at point B ( X c - 10 -3 ) after 9 Gyr, after which nuclear 
energy generation gradually moves out to a thick shell surrounding the isothermal helium core. Be¬ 
tween B and C the core slowly grows in mass and contracts, while the envelope expands in response 
and the burning shell gradually becomes thinner in mass. By point C the helium core has become 
degenerate. At the same time the envelope has cooled and become largely convective, and the star 
finds itself at the base of the red giant branch (RGB), close to the Hayashi line. The star remains 
in thermal equilibrium throughout this evolution and phase B-C lasts about 2 Gyr for this 1 M© star. 
This long-lived phase corresponds to the well-populated subgiant branch in the H-R diagrams of old 
star clusters. 

Stars with masses in the mass range 1.1 - 1.5 M© show a very similar behaviour after the main 
sequence, the only difference being the small convective core they develop during core H-buming. 
This leads to a ‘hook’ in the evolution track at central H exhaustion (see Sec. 8.3). The subsequent 
evolution during H-shell burning is similar, the core remaining in TE until it becomes degenerate on 
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Figure 9.5. Evolution of a 1 M e star of ini¬ 
tial composition X = 0.7, Z = 0.02. The 
top panel (a) shows the internal structure as 
a function of mass coordinate m. Gray areas 
are convective, lighter-gray areas are semi- 
convective. The red hatched regions show 
areas of nuclear energy generation: e nuc > 
5L/M (dark red) and e nuc > L/M (light 
red). The letters A... J indicate correspond¬ 
ing points in the evolution track in the H-R 
diagram, plotted in the bottom panel (b). See 
text for details. 


the RGB and a correspondingly slow evolution across the subgiant branch. Stars with 1.5 < M/M 0 < 
2 do exhibit a small Hertzsprung gap as they reach the Schonberg-Chandrasekhar limit before their 
cores become degenerate. After a period of slow, thick shell burning on the subgiant branch they 
undergo a phase of rapid, thermal-timescale expansion until they reach the giant branch. In this case 
the gap in T c a to be bridged is narrow because the main sequence is already relatively close in effective 
temperature to the Hayashi line. 

Regardless of these differences between stars of different mass during the early shell-H burning 
phase, all stars with M < 2 M 0 have in common that their helium cores become degenerate before 
the central temperature is high enough for helium ignition, and they settle into TE on the red giant 
branch. 
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9.2.3 The red giant branch in low-mass stars 

The evolution of low-mass stars along the red giant branch is very similar and almost independent 
of the mass of the star. The reason for this similarity is that by the time the helium core has become 
degenerate, a very strong density contrast has developed between the core and the envelope. The 
envelope is so extended that it exerts very little weight on the compact core, while there is a very 
large pressure gradient between core and envelope. The pressure at the bottom of the envelope (see 
eq. 8.14) is very small compared to the pressure at the edge of the core and in the hydrogen-burning 
shell separating core and envelope. Therefore the stellar structure depends almost entirely on the 
properties of the helium core. Since the core is degenerate, its structure is independent of its thermal 
properties (temperature) and only depends on its mass. Therefore the structure of a low-mass red 
giant is essentially a function of its core mass. 

As a result there is a very tight relation between the helium core mass and the luminosity of a red 
giant, which is entirely due to the hydrogen shell-burning source. This core-mass luminosity relation 
is very steep for small core masses, M c < 0.5 M© and can be approximately described by a power law 

L«2.3xl0 5 L©(—) (9.2) 

\ M eJ 

Note that the luminosity of a low-mass red giant is independent of its total mass. Therefore the 
evolution of all stars with M < 2 M© converges, and also the central density and temperature start 
following the same evolution track (e.g. see Fig. 7.4). In the H-R diagram the star is located along 
the Hayashi line appropriate for its mass M, meaning that higher-mass red giants have slightly higher 
T e ff at the same luminosity. (This means there is also a core-mass radius relation, but it is less tight 
than the M c -L relation and depends slightly on the total mass.) 

As the H-buming shell adds mass to the degenerate helium core, the core slowly contracts and the 
radius and luminosity increase. The higher luminosity means the H-shell must bum at a higher rate, 
leading to faster core-mass growth. The evolution along the RGB thus speeds up as the luminosity 
increases (see Fig. 9.5). The density contrast between core and envelope increases and the mass within 
the burning shell decreases, to « 0.001 M© near the tip of the RGB. Since less mass is contained 
in the burning shell while the luminosity increases, the energy generation rate per unit mass e nuc 
increases strongly, which means the temperature within the burning shell also increases. With it, the 
temperature in the degenerate helium core increases. When the tip of the RGB is reached (at point 
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Figure 9.6. Evolution track of a 0.8 M 0 star of rather 
low metallicity, Z = 0.004. The inset shows the tem¬ 
porary decrease of luminosity when the H-burning 
shell crosses the hydrogen discontinuity left by the 
first dredge-up (corresponding to point E in Fig. 9.5). 
The open circle indicates where first dredge-up oc¬ 
curs. Figure from Salaris & Cassisi. 
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F in Fig. 9.5) at L « 2000 L Q and a core mass of ~ 0.45 M 0 , the temperature in the degenerate core 
has reached a value close to 10 8 K and helium is ignited. This is an unstable process due to the strong 
degeneracy, and leads to a thermonuclear runaway known as the helium flash (see Sec. 9.3.2). 

First dredge-up and the luminosity bump 

When the convective envelope reaches its deepest extent at point D in Fig. 9.5, it has penetrated into 
layers that were processed by H-buming during the main sequence, and have been partly processed 
by the CN-cycle. Up to point D the surface He abundance increases and the H abundance decreases, 
but more noticeably the C/N ratio decreases by a large factor. This is called the first dredge-up phase 
(later dredge-ups occur after He burning). 

Some time later, at point E in Fig. 9.5 the H-buming shell has eaten its way out to the discontinuity 
left by the convective envelope at its deepest extent. The shell suddenly finds itself in an environment 
with a higher H abundance (and a lower mean molecular weight). As a consequence it starts burning 
at a slightly lower rate, leading to a slight decrease in luminosity (see Fig. 9.6). The resulting loop 
(the star crosses this luminosity range three times) results in a larger number of stars in this luminosity 
range in a stellar population. This ‘bump’ in the luminosity function has been observed in many old 
star clusters. 


Another process that becomes important in low-mass red giants is mass loss. As the stellar luminosity 
and radius increase as a star evolves along the giant branch, the envelope becomes loosely bound and 
it is relatively easy for the large photon flux to remove mass from the stellar surface. The process 
driving mass loss in red giants is not well understood. When calculating the effect of mass loss in 
evolution models an empirical formula due to Reimers is often used: 

M = -4xlO- 13 77 ^£^M 0 /yr (9.3) 

where rj is a parameter of order unity. It should be noted that this relation is based on observations 
of only a handful of stars with well-determined stellar parameters. A value of rj ~ 0.25 - 0.5 is often 
used because it gives the right amount of mass loss on the RGB to explain the morphology in the H-R 
diagram of stars in the subsequent helium-burning phase, on the horizontal branch. The 1 M 0 star of 
our example loses about 0.3 M 0 of its envelope mass by the time it reaches the tip of the giant branch. 

9.3 The helium burning phase 

As the temperature in the helium core approaches 10 s K, the 3a reaction starts to produce energy at 
a significant rate. This is the onset of the helium burning phase of evolution. Unlike for hydrogen 
burning, the reactions involved in helium burning (see Sect. 5.4.4) are the same for all stellar masses. 
However, the conditions in the core at the ignition of helium are very different in low-mass stars 
(which have degenerate cores) from stars of higher mass (with non-degenerate cores). Therefore 
these cases will be discussed separately. 

9.3.1 Helium burning in intermediate-mass stars 

We again take the 5 M 0 star depicted in Figs. 9.2-9.4 as a typical example of an intermediate-mass 
star. The ignition of helium takes place at point E in these figures. Since the core is non-degenerate at 
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this point (p c ~ 10 4 g/cm 3 , Fig. 9.1), nuclear burning is thermally stable and helium ignition proceeds 
quietly. Owing to the high temperature sensitivity of the He-buming reactions, energy production 
is highly concentrated towards the centre which gives rise to a convective core. The mass of the 
convective core is 0.2 M© initially and grows with time (unlike was the case for hydrogen burning). 

Initially, the dominant reaction is the 3 a reaction which converts 4 He into 12 C inside the convec¬ 
tive core. As the 12 C abundance builds up, the 12 C +a reaction gradually takes over, so that 16 0 is also 
produced at a rate that increases with time (see Fig. 9.4b and compare to Fig. 5.5). When the central 
He abundance X He < 0.2 the mass fraction of 12 C starts decreasing as a result of the diminishing 
3 a rate (which is proportional to X 3 e ). The final 12 C/ 16 0 ratio is about 0.3, decreasing somewhat 
with stellar mass. This is related to the fact that in more massive stars the central temperature during 
He burning is larger. Note that the final 12 C/ 16 0 ratio depends on the uncertain rate of the l2 C(a, y) 
reaction, and the values given here are for the rate that is currently thought to be most likely. 

The duration of the central helium burning phase in our 5 M© star (E-H) is about 22Myr, i.e. 
approximately 0.27 X tms- This seems surprisingly long given that the energy gain per gram of He 
burning is only 10 % of that of H burning, while the luminosity of the star is (on average) somewhat 
larger than during the main sequence. The reason can be discerned from Fig. 9.4a: most of the 
luminosity during helium burning still comes from the H-burning shell surrounding the core, although 
the luminosity contribution of He burning (T He ) increases with time and becomes comparable towards 
the end of this phase. 

We can understand the behaviour of T Hc by considering that the properties of the helium core 
essentially depend only on the core mass M c and are hardly affected by the surrounding envelope. 
Because the envelope is very extended the pressure it exerts on the core (eq. 8.14) is negligible com¬ 
pared to the pressure inside the dense helium core. In fact is a steep function of M c , analogous to 
the main-sequence M-L relation - indeed, if the envelope were stripped away, the bare helium core 
would lie on a helium main sequence? As a result of H-shell burning, M c grows with time during the 
He-burning phase and T Hc increases accordingly. Another consequence is that in models computed 
with convective overshooting L Hc is larger on account of the larger core mass left after the main se¬ 
quence (see Sect. 8.3.4). Therefore the duration of the He burning phase (i.e. the appropriate nuclear 
timescale, r nuc M c /Ln e ) is shorter in models with overshooting. A 5 M© star of the same composi¬ 
tion computed with overshooting has a main-sequence lifetime tms = 100 Myr and a helium-burning 
lifetime of 16 Myr. 

During helium burning intermediate-mass stars describe a loop in the H-R diagram (E-H in 
Fig. 9.2). After He ignition at the tip of the giant branch, the envelope contracts (on the nuclear 
timescale for helium burning) and the stellar radius decreases. Initially the luminosity also decreases 
while the envelope is mostly convective (E-F) and the star is forced to move along its Hayashi line. 
When most of the envelope has become radiative at point F, the star leaves the red giant branch and 
the effective temperature increases. This is the start of a so-called blue loop, the hottest point of which 
is reached at G when X He ~ 0.3. This also corresponds to a minimum in the stellar radius, after which 
the envelope starts expanding and the star again approaches the giant branch when X He ~ 0.05. By 
the end of core helium burning (H) the star is back on the Hayashi line, very close to its starting point 
(E). If we consider stars of different masses, the blue extension of the loops in the HRD increases (the 
loops extend to larger T e s values) for increasing mass, up to M « 12 M©. (The behaviour of stars of 
larger masses can be more complicated, one of the reasons being strong mass loss, and we defer a 
discussion of this until Chapter 11.) On the other hand, for M < 4 M© the loops always stay close to 
the red giant branch and do not become ’blue’. 

The occurrence of blue loops is another example of a well-established result of detailed stellar 

2 The mass-luminosity relation for such helium main-sequence stars can be approximately described by the homology 
relation (6.32) if the appropriate value of u is used. 
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Figure 9.7. The helium flash. Evolution with time of the surface luminosity ( L s ), the He-burning luminosity 
{L'ia) and the H-buming luminosity (Lh) during the onset of He burning at the tip of the RGB in a low-mass 
star. Time t = 0 corresponds to the start of the main helium flash. Figure from Salaris & Cassisi. 


evolution calculations, that is difficult to explain in terms of basic physics. The detailed models 
indicate that the occurrence and extension of blue loops depends quite sensitively on a number of 
factors: the chemical composition (mainly Z), the mass of the helium core relative to the envelope, 
and the shape of the hydrogen abundance profile above the core. It therefore also depends on whether 
convective overshooting was assumed to take place during the main sequence: this produces a larger 
core mass, which in turn has the effect of decreasing the blue-ward extension of the loops while 
increasing their luminosity. 

The blue loops are important because they correspond to a slow, nuclear timescale phase of evo¬ 
lution. One therefore expects the corresponding region of the H-R diagram to be well populated. 
More precisely, since intermediate-mass stars spend part of their He-buming phase as red giants and 
part of it in a blue loop, one expects such stars to fill a wedge-shaped region in the HRD. Indeed one 
finds many stars in the corresponding region, both in the solar neighbourhood (Fig. 1.1, although this 
is dominated by low-mass stars) and in open clusters with ages less than ~ 1 Gyr. The dependence 
of the loops on overshooting also makes observational tests of overshooting using He-buming stars 
possible. Another significant aspect of blue loops is that they are necessary for explaining Cepheid 
variables (see Sect. 9.3.3), which are important extragalactic distance indicators. 

9.3.2 Helium burning in low-mass stars 

In low-mass stars (with M <2 M 0 ) the helium burning phase differs from more massive stars in two 
important aspects: (1) helium ignition occurs under degenerate conditions, giving rise to a helium 
flash, and (2) all low-mass stars start helium burning with essentially the same core mass M c ~ 
0.45 M q (Sect. 9.2.3). The luminosity of low-mass He-buming stars is therefore almost independent 
of their mass, giving rise to a horizontal branch in the HRD. 

The helium flash 

We again take a star of 1 M e as a typical example of all low-mass stars. Helium ignition occurs 
when T c ~ 10 8 K and p c ~ 10 6 g/cm 3 , so the helium core is strongly degenerate (see Fig. 9.1). 
We have seen in Sect. 6.5.2 that helium burning under these conditions is thermally unstable: the 
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energy generated by the 3a reaction causes a temperature increase, rather than a decrease, and helium 
ignition thus initiates a thermonuclear runaway. The reason is that the degenerate pressure is basically 
independent of T, so that the energy released by fusion does not increase the pressure and therefore 
leads to negligible expansion and negligible work done. All nuclear energy released therefore goes 
into raising the internal energy. Since the internal energy of the degenerate electrons is a function 
of p and hence remains almost unchanged, it is the internal energy of the non-degenerate ions that 
increases and thus raises the temperature. As a result, the evolution is vertically upward in the p c -T c 
diagram. 3 

The thermonuclear runaway leads to an enormous overproduction of energy: at maximum, the 
local luminosity in the helium core is l ~ 10 10 L e - similar to a small galaxy! However, this only lasts 
for a few seconds. All the nuclear energy released is absorbed by expansion of the non-degenerate 
layers surrounding the core, so none of this luminosity reaches the surface. Since the temperature 
increases at almost constant density, degeneracy is eventually lifted when T « 3 x 10 s K. Further 
energy release increases the pressure when the gas starts behaving like an ideal gas and thus causes 
expansion and cooling. This results in a decrease of the energy generation rate until it balances the 
energy loss rate and the core settles in thermal equilibrium at T c w 10 s K and p c ~ 2 x l 0 4 g/cm 3 (see 
Fig. 9.1). Further nuclear burning is thermally stable. 

Detailed numerical calculations of the helium flash indicate that this sequence of events indeed 
takes place, but helium is not ignited in the centre but in a spherical shell at m « 0.1 M Q where T has 
a maximum. This off-centre temperature maximum is due to neutrino losses during the preceding red 
giant phase. At high density and temperature, neutrinos are released by a number of weak interaction 
processes (see K&W Chapter 18.6). Since neutrinos thus created escape without interacting with the 
stellar gas, this leads to effective cooling of the central region of the degenerate helium core. The 
mass coordinate at which T max occurs (and where helium ignites) decreases somewhat with stellar 
mass. 

The high local luminosity causes almost the entire region between the the ignition point (at m ~ 
0.1 M q ) up to the bottom of the H-buming shell (at 0.45 M©) to become convective. The energy 
released in the He flash is thus transported efficiently to the edge of the core, where it is absorbed 
by expansion of the non-degenerate layers. Convection also mixes the product of the He flash ( 12 C 
produced in the 3a reaction) throughout the core. About 3 % of the helium in the core is converted 
into carbon during the flash. Because the convective shell containing this carbon never overlaps 
with the convective envelope surrounding the H-buming shell, this carbon does not reach the surface. 
(However, this may be different at very low metallicity.) 

After the He flash, the whole core expands somewhat but remains partially degenerate. In detailed 
models a serious of smaller flashes follows the main He flash (see Fig. 9.7) during « 1.5 Myr, before 
degeneracy in the centre is completely lifted and further He burning proceeds stably in a convective 
core, as for intermediate-mass stars. 

The horizontal branch 

In our 1 Mq example star, the helium flash occurs at point F in Fig. 9.5. Evolution through the helium 
flash was not calculated for the model shown in this figure. Instead, the evolution of the star is resumed 
at point G when the helium core has become non-degenerate and has settled into TE with stable He 
burning in the centre and H-shell burning around the core. (Models constructed in this way turn out 
to be very similar to models that are computed all the way through the He flash, such as shown in 
Fig. 9.7.) At this stage the luminosity and radius of the star have decreased by more than an order of 

3 This part of the evolution is skipped in the 1 M 0 model shown in Fig. 9.1, which is why a gap appears in the evolution 
track. The evolution during the He flash is shown schematically as a dashed line for the 1 M 0 model in Fig. 7.4. 
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Figure 9.8. Location of the zero-age 
horizontal branch (think gray line) for a 
metallicity Z = 0.001 typical of glob¬ 
ular clusters. These models have the 
same core mass (0.489 M 0 ) hut varying 
total (i.e. envelope) mass, which deter¬ 
mines their position in the H-R diagram. 
Evolution tracks during the HB for sev¬ 
eral total mass values are shown as thin 
solid lines. Figure from Maeder. 


magnitude from the situation just before the He flash. Here we again see the mirror principle at work: 
in this case the core has expanded (from a degenerate to a non-degenerate state) and the envelope has 
simultaneously contracted, with the H-buming shell acting as a ‘mirror’. 

In the 1 M q star of solar composition shown in Fig. 9.5, helium burning occurs between G and 
H. The position of the star in the H-R diagram does not change very much during this period, always 
staying close (but somewhat to the left of) the red giant branch. The luminosity is « 50 L© for most 
of the time; this value is determined mainly by the core mass. Since the core mass at the start of 
helium burning is ~ 0.45 M© for all low-mass stars, independent of stellar mass, the luminosity at 
which He burning occurs is also almost independent of mass. If we consider He-buming stars of 
a given composition (e.g. in a star cluster), only the envelope mass may vary from star to star. At 
solar metallicity, all such stars occupy about the same position in the HRD. This gives rise to a so- 
called red clump in observed colour-magnitude diagrams of low-mass stellar populations (visible for 
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Figure 9.9. Evolution with time of the lu¬ 
minosities and central abundances in a 1M 0 
star during the late part of the red giant branch 
and during helium burning. Letters D... H 
correspond to the same evolution phases as 
in Fig. 9.5. 


age (10 9 yr) 
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instance in Fig. 1.1). However, the radius and effective temperature of He-buming stars depends on 
their envelope mass. Stars with a small envelope mass (either because of a smaller initial mass, or 
because they suffered a larger amount of mass loss on the RGB) can be substantially hotter than the 
one shown in Fig. 9.5. Furthermore, at low metallicity the critical envelope mass, below which He- 
buming stars become small and hot, is larger. Stars with different amounts of mass remaining in their 
envelopes can then form a horizontal branch in the HRD (Fig. 9.8). Horizontal branches are found in 
old stellar populations, especially in globular clusters of low metallicity (an example is the globular 
cluster M3 shown in Fig. 1.2). The observed distribution of stars along the HB varies greatly from 
cluster to cluster, and the origin of these different HB morphologies is not fully understood. 

The duration of the core helium burning phase is about 120 Myr, again independent of stellar 
mass. While this is longer than in intermediate-mass stars, it is a much shorter fraction of the main- 
sequence lifetime because of the much higher luminosity of the He-buming phase. The evolution of 
the stellar structure during helium burning is qualitatively similar to that of intermediate-mass stars; 
see Figs. 9.5a and 9.9. The most striking differences are: 

• The contribution of He-buming to the stellar luminosity is larger, especially towards the end of 
the phase. This is due to the relatively small envelope mass. 

• The development of a substantial semi-convective region on top of the convective core. This 
is related to a difference in opacity between the C-rich convective core and the He-rich zone 
surrounding it, and gives rise to partial (non-homogeneous) mixing in this region. 

• The occurrence of ‘breathing pulses’, giving rise to the sudden jumps in the central composition 
and in the luminosity. Whether these are real or simply a numerical artifact of one-dimensional 
stellar models is not clear. 4 

9.3.3 Pulsational instability during helium burning 

Cepheid variables are luminous pulsating stars with periods between about 2 and 100 days. It turns 
out that there is a well-defined correlation between the pulsation period and the luminosity of these 
stars, first discovered for Cepheids in the Small Magellanic Cloud. A modem version of this empirical 
relation is shown in Fig. 9.11. Their importance for astronomy lies in the fact that the period can be 
easily determined, even for stars in other galaxies, and thus provides an estimate of the absolute 
luminosity of such a star, making Cepheids important standard candles for the extragalactic distance 
scale. 

Cepheids lie along a pulsational instability strip in the H-R diagram (see Fig. 9.10). During 
the evolution of an intermediate-mass star, this instability strip is crossed up to three times. The 
first crossing occurs during H-shell burning (C-D in Fig. 9.2) but this is such a rapid phase that the 
probability of catching a star in this phase is very small. In stars with sufficiently extended blue loops, 
another two crossings occur (F-G and G-H) during a much slower evolution phase. Cepheids must 
thus be helium-burning stars undergoing a blue loop. 

Since pulsation is a dynamical phenomenon, the pulsation period is closely related to the dy¬ 
namical timescale (eq. 2.17). Therefore the pulsation period n is related the mean density: to first 
approximation once can write n oc p~ 1 1 1 oc M -1 / 2 /? 3 / 2 . Each passage of the instability strip yields a 
fairly well-defined radius and luminosity. Passage at a larger L corresponds to a larger R and therefore 
to a larger n, because the variation in mass is smaller than that in radius and enters the relation with 

4 For details about the latter two effects, see either Salaris & Cassisi or John Lattanzio’s tutorial at 
http://www.maths.monash.edu.au/johnl/StellarEvolnDemo/. 
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Figure 9.10. Occurrence of various classes 
of pulsating stars in the H-R diagram, over¬ 
laid on stellar evolution tracks (solid lines). 
Cepheid variables are indicated with ‘Ceph’, 
they lie within the pulsational instability strip 
in the HRD (long-dashed lines). Their equiv¬ 
alents are the RR Lyrae variables among 
HB stars (the horizontal branch is shown as 
a dash-dotted line), and the 6 Scuti stars 
(6 Set) among main-sequence stars. Pulsa¬ 
tional instability is also found among lumi¬ 
nous red giants (Mira variables), among mas¬ 
sive main-sequence stars - ySCep variables 
and slowly pulsating B (SPB) stars, among 
extreme HB stars known as subdwarf B stars 
(sdBV) and among white dwarfs. Figure 
from Christensen-Dalsgaard (2004). 


a smaller power. This provides a qualitative explanation of the period-luminosity relation. The min¬ 
imum observed period should correspond to the the lowest-mass star undergoing a blue loop. Also 
the number of Cepheids as a function of period must correspond to the time it takes for a star of the 
corresponding mass to cross the instability strip. Thus Cepheids provide a potential test of stellar 
evolution models. 


Suggestions for further reading 

The contents of this chapter are also covered by Chapters 25.3.2 and 26.1-26.5 of Maeder, while 
stellar pulsations and Cepheids are treated in detail in Chapter 15. See also Kippenhahn & Weigert, 
Chapters 31 and 32. 
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Figure 9.11. The period-luminosity relation for classical Cepheids in the Large Magellanic Cloud. Luminosity 
is expressed as absolute magnitude in the B band (left) and in the V band. Figure from Sandage et al. (2004, 
A&A 424, 43). 

Exercises 


9.1 Conceptual questions 

(a) Why does the luminosity of a star increase on the main sequence? Why do low-mass stars, like 
the Sun, expand less during the main sequence than higher-mass stars? 

(b) Explain what happens during the ‘hook’ at the end of the main sequence of stars more massive 
than the Sun. 

(c) What is convective overshooting ? Think of at least three effects of overshooting on the evolution 
of a star. 

(d) Explain the existence of a Hertzsprung gap in the HRD for high-mass stars. Why is there no 
Hertzsprung gap for low-mass stars? 

(e) What do we mean by the mirror principle ? 

(f) Why does the envelope become convective on the red giant branch? What is the link with the 
Hayashi line ? 

9.2 Evolution of the abundance profiles 

(a) Use Fig. 9.4 to sketch the profiles of the hydrogen and helium abundances as a function of the 
mass coordinate in a 5 M 0 star, at the ages corresponding to points C, E, G and H. Try to be as 
quantitative as possible, using the information provided in the figure. 

(b) Do the same for a 1 M 0 star, using Figs. 9.5 and 9.9, at points B, D, F and H. 

(c) The abundances plotted in Figs. 9.4 and 9.9 are central abundances. What happens to the abun¬ 
dances at the surface? 

9.3 Red giant branch stars 

(a) Calculate the total energy of the Sun assuming that the density is constant, i.e. using the equation 
for potential energy £ gr = GM 2 /R. In later phases, stars like the Sun become red giants, with 
R « 100/? Q . What would be the total energy, if the giant had constant density. Assume that the 
mass did not change either. Is there something wrong? If so, why is it? 
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(b) What really happens is that red giants have a dense, degenerate, pure helium cores which grow to 
~ 0.45M o at the end of the red giant branch (RGB). What is the maximum radius the core can 
have for the total energy to be smaller than the energy of the Sun? (N.B. Ignore the envelope - 
why are you allowed to do this?) 

(c) For completely degenerate stars, one has 


R = 2.6 


10 


(9.4) 


where p e is the molecular weight per electron and p e = 2 for pure helium. Is the radius one finds 
from this equation consistent with upper limit derived in (b)? 


9.4 Core mass-luminosity relation for RGB stars 

Low-mass stars on the RGB obey a core mass-luminosity relation, which is approximately given by 
eq. (9.2). The luminosity is provided by hydrogen shell burning. 


(a) Derive relation between luminosity L and the rate at which the core grows dMJdt. Use the energy 
released per gram in hydrogen shell burning. 

(b) Derive how the core mass evolves in time, i.e, M c = M c (t). 

(c) Assume that a star arrives to the RGB when its core mass is 15% of the total mass, and that it 
leaves the RGB when the core mass is 0.45 M 0 . Calculate the total time a 1 M 0 star spends on the 
RGB and do the same for a 2 M e star. Compare these to the main sequence (MS) lifetimes of these 
stars. 

(d) What happens when the core mass reaches 0.45 M 0 ? Describe the following evolution of the star 
(both its interior and the corresponding evolution in the HRD). 

(e) What is the difference in evolution with stars more massive than 2 M 0 ? 


9.5 Jump in composition 


Consider a star with the following distribution of hydrogen: 


X(m) = 


0.1 for m < m c 
0.7 for m < m c 


(9.5) 


(a) In this star a discontinuous jump in the composition profile occurs at m = m c . What could have 
caused such a chemical profile? Explain why P and T must be continuous functions. 

(b) Calculate the jump in density Ap/p. 

(c) Also calculate the jump in opacity, A k/k, if the opacity is given as: 

- Kramers: /q,/ ~ Z( 1 + X)pT~ 3 5 

- Electron scattering: K e = 0.2(1 + X) 


137 



